Ionization wave propagation on a micro cavity plasma array 
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Abstract 

Microcavity plasma arrays of inverse pyramidal cavities have been fabricated in p-Si wafers. 
Each cavity acts as a microscopic dielectric barrier discharge. Operated at atmospheric pressure in 
argon and excited with high voltage at about 10 kHz, each cavity develops a localized microplasma. 
Experiments have shown a strong interaction of individual cavities, leading to the propagation of 
wavedike optical emission structures along the surface of the array. This phenomenon is numerically 
investigated using computer simulation. The observed ionization wave propagates with a speed of 
about 5 km/s, which agrees well the experimental findings. It is found that the wave propagation 
is due to sequential contributions of a drift of electrons followed by drift of ions between cavities 
seeded by photoemission of electrons by the plasma in adjacent cavities. 



Atmospheric pressure microplasmas have gained increasing attention since the early work 
on microhollow cathode discharges by Schoenbach et al. [1]. Although they are highly 
collisional, microplasmas show strong non-equilibrium behavior [2]. Electrons are hot (up 
to a few eV) while ions and and the neutral gas are cold (about room temperature in most 
cases, though temperatures of many hundreds Centigrade occur at high power loading). 
Microplasmas are characterized by high electron densities (up to a few 10 17 cm -3 ) and high 
power densities (up to some 100 kW/cm 3 ). Due to these unique properties microplasmas 
have found widespread use in technological and biomedical applications [3-6] , such as surface 
treatment, sterilization, and lighting. At the same time microplasmas continue to show 
unexpected phenomena and are therefore scientifically interesting. For example, arrays 
of microplasmas having thousands of individual microdischarges have been fabricated as 
efficient lighting sources [7]. When excited by AC voltages, the initiation of optical emission 
from such arrays has shown wave-like phenomena, propagating across the array at speeds 
of a few kilometers per second. In this paper, we discuss results from experimental and 
computational investigations of these wave-like propagation of optical emission. 

The microcavity array we investigated, developed by Eden et al. [7], has individual 
inverse pyramidal cavities fabricated on p-Si wafers (see Fig. 1), combined with a nickel 
grid as a counter electrode and dielectric coatings. Each cavity - in essence a microscopic 
dielectric barrier discharge (DBD) - has a base opening of 50 x 50 /im and a depth of 35 
fim. The separation of the individual cavities is 50 /im. The entire array consists of 50 by 
50 microplasmas. When operated with a sawtooth voltage of about 10 kHz in atmospheric 
pressure argon, the array emits a bright glow in the visible wavelength range, which visually 
appears homogeneous over the entire array. However, spatially and temporally resolved 
emission spectroscopy shows that wave-like structures (see Fig. 2) propagate across the 
array, indicating that the individual microdischarges strongly interact with each other [8, 9]. 

The related computer simulation was performed using the modeling platform nonPDP- 
SIM, described in detail in [10-12]. Here we briefly discuss the relevant physical equa- 
tions. Poissons equation for the electrostatic potential is self-consistently coupled with 
drift-diffusion equations for the charged species and the surface charge balance equation. 
The set of equations is simultaneously integrated in time using an implicit Newton iteration 
technique. This integration step is followed by an implicit update of the electron temper- 
ature by solving the electron energy equation. To capture the non-Maxwellian behavior of 
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the electrons, the electron transport coefficients and rate coefficients are obtained by solving 
the zero-dimensional Boltzmann's equation for the electron energy distribution. A Monte 
Carlo simulation is used to track the trajectories of sheath accelerated secondary electrons. 
The transport of photons is treated by means Green's function method. The discharge is 
sustained in argon at atmospheric pressure. The species in the model are electrons, Ar(3s), 
Ar(4s), Ar(4p), Ar + , Ar£ and Ar^~. The photon transport we tracked in the model is emis- 
sion dimmer radiation Ar^. In the absence of impurities and large densities of excited states, 
this emission is optically thin. The reaction mechanism is summarized in Ref. [13]. 

The simulation addresses a sub-system of the micro cavity array over a time interval of 
50 ns. The geometry is two-dimensional and consists of three neighboring cavities, which 
appear in the model as slits. The entire computational domain is 1 mm tall by 1.5 mm 
wide, with boundaries of which are electrically grounded. The conductivity of the wafer 
in contact with the boundary is 0.17 ficm. The secondary electron emission coefficient by 
ions is 0.15 on all surfaces. A photoelectron emission coefficient of 0.05 is assumed on the 
dielectric surface. The step function DC potential of —500 V is applied to the embedded 
nickel electrode grid. 

The initial plasma density is a Gaussian shaped spot of plasma 10 fim wide and 5 x 10 10 
cm -3 in magnitude at the left vertex of the left cavity. The resulting electron density as 
a function of time is shown in Fig. 3. Ions accelerated by the negative potential of the 
nickel grid drift towards the vertex of the cavity and initiate secondary electron emission. 
The secondary electrons are accelerated into the cavity and generate ions and electrons. The 
ionization rate peaks at the vertices at about 5 x 10 24 cm -3 s _1 with an electron temperature 
of 8 eV, and then spreads on the top dielectric, similar to a conventional DBD. Meanwhile, 
photons produced by relaxation of Ar^j lead to photoelectron emission on the dielectric 
surfaces, mainly at the right vertex of the neighboring cavity. At t = 25 ns a sufficiently 
large electron density of about 10 15 cm -3 is reached in the left cavity, that the potential is 
shielded from the interior of the left cavity and the discharge turns from a Townsend mode 
to a glow mode. The quasineutral region fills the left cavity and spreads over its edges, 
charging the top dielectric. 5 ns after the left cavity is ignited, an electron density larger 
than 10 10 cm~ 3 is produced in the middle cavity initiated by avalanching of photo-generated 
electrons, and another 5 ns later in the right cavity. Both cavities undergo a transition from a 
Townsend mode to glow mode 10 ns after the initial electron density is generated (Fig. 3d,e). 
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In all three cavities a quasineutral region develops. These individual plasmas connect after 
their expansion out of the cavities, as shown in Fig. 3f. The electron density peaks in the 
head of the avalanche with n e = 1.5 x 10 16 cm -3 . Within the cavities the electron density 
exceeds 10 15 cm -3 leading to a Debye length Ad < 1 jt/m. During the ignition, the electron 
temperature reaches 5 eV within the cavities, but decays rapidly. The observed densities 
agree with values expected for this kind of discharge [3]. 

The electric potential is plotted at t = 40 ns in Fig. 4. It shows the three phases of the 
discharge development. The potential in the right cavity is essentially the vacuum potential 
- there is little space charge to influence the electrical field. The left cavity in contrast 
contains a quasineutral plasma with the plasma being nearly equipotential. The majority of 
the potential is dropped over the dielectrics which covers the anode and the cathode, as one 
expects for a dielectric barrier discharge. The electron avalanche within the center cavity 
produces an equipotential channel which merges with that of the left cavity and indicates 
the transition from a Townsend mode to a glow mode discharge. 

Taking the time between the ignition of the left and the right cavity one finds an effective 
ionization propagation speed of about 5 km/s which agrees very well with the measured 
propagation speed of 5 to 10 km/s. In the absence of photoelectron emission, the propaga- 
tion speed is about 1 km/s. Although not all effects are addressed with our simulation, the 
ionization wave is well described. Other effects of importance are surface charges accumu- 
lated in a previous half cycle and a second ionization wave after an further increase of the 
voltage. The effect of surface charges accumulated at a previous half cycle with an inverse 
voltage results in an enhancement of the electrical field an reduces the breakdown voltage. 
Both are a matter of ongoing research. To summarize, the observed ionization wave is driven 
by photoelectrons emitted at neighbored cavities. These electron source causes the ignition 
of a single cavity and provides new photons to sustain the ionization wave propagation. 

The authors gratefully acknowledge financial support by the Deutsche Forschungsgemein- 
schaft in the frame of Research Group 1123 Physics of Microplasmas as well as the Ruhr 
University Research School. 
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FIG. 2: Experimental result from phase resolved optical emission spectroscopy of the microplasma 
array as registered with an ICCD camera (false colors) for a gate interval of 100 ns and spectrally 
integrated, (p = 10 5 Pa Ar, U PP = 830 V, / ac = 20 kHz.) Image taken at V{t) = -400 V applied 
to Ni-grid. 
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FIG. 3: Temporal evolution of the electron density distribution in units of cm" 3 on logarithmic 
scale. Figures (a-f) shows the ignition of the three cavities starting at 25 ns in time steps of At 
= 5 ns. The width and separation of a single cavity is 50 //m, the overall depth is 45 /jiii. The 
simulation is performed for Ar as a working gas at atmospheric pressure, a driving DC voltage of 
—500 V applied to the embedded Ni-grid. 
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FIG. 4: Simulated potential at t = 80 ns. Left hand cavity is ignited whereas right hand cavity 
shows vacuum potential. Centered cavity is in the stage of transition from Townsend mode to glow 
mode. 



